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Abstract 

Dihadron azimuthal correlation measurements have revealed striking modifications of the jets 
by the dense medium created in heavy-ion collisions at RHIC. One important question is to what 
extent the modification can be attributed to cold nuclear matter effects. In this analysis, we carried 
out a detailed mapping of the correlation patterns using high-statistics RUNS d+Au minimum bias 
data. A striking scaling behavior of the jet pair yields is observed at low and intermediate pT- The 
jet pair yields are found to be enhanced relative to A'^coii scaled p-|-p jet pair yields. The nuclear 
modification factor for jet pair yields, JdAu, seems to scale with p^'^= p9p + p\' (scaler sum), and 
shows a characteristic Cronin-like enhancement at <5-7GeV/c. Interestingly, the level of yield 
modifications is similar between the near- and away-side pairs, and the jet shapes are not modified 
relative to p+p collisions. The pedestal yield under the jet peak is studied in p+p and d+Au 
collisions. The pedestal yield in p+p collisions is found to be larger than PYTHIA calculations. 
In d+Au collisions, it is found to exceed a simple sum of one p+p jet event and Ncoii — 1 minimum 
bias p+p events. The possible interpretation of these results and their implications for Au+Au 
measurements are discussed. 

PACS numbers: 25.75.Dw 13.87.Fh 
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I. INTRODUCTION 



It is generally believed that a strongly interacting Quark Gluon Plasma (sQGP) is created 
in central Au+Au collisions at the relativistic heavy ion collider (RHIC). The two most 
important evidences for sQGP are the observation of a large elliptic flow and a strong 



suppression of high jets (jet quenching 
measurements of single particle production 



various correlation measurements 



flQ. 




hese results were initially obtained from 
and were subsequently confirmed by 



The information implied by correlation measurements, however, are much richer than by 
single particle measurements. The complicated correlation patterns can be better under- 
stood by studying them separately in a high px region and a low region. The high p^ 
region (> 5 GeV/c) is dominated by a suppressed but essentially p+p like jet fragmentation 
component, which can be interpreted as the combined result of jet quenching and surface 
bias [9]. The low px region is characterized by a highly non-trivial modification that de- 
pends on the pt, A(/) and Arj. This latter region seems to be driven by a detailed balance 
between the jet and the medium: the "jet" signal exhibits the famous "ridge" and "cone" 
like structures in Arj and A0 j^, Q]; the bulk medium manifests a quark number scaling 
of elliptic fiow 10[ and an enhanced baryon/meson ratio Current efforts are focused 
on obtai ning a quantitative understanding of the mutual influences between the jet and the 



medium 
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13| . Both the modiflcation of jets by the flowing medium and the response of 



the medium to quenched jets need to be taken in to account. 

Study of single hadron production and dihadron correlation in p+A or p+A-like collision 
such as d+Au provides important baselines for understanding the results obtained in Au+Au 
collisions. Previous measurements established the dominance of the flnal state effects for 
the modiflcations of single particle production in Au+Au collisions [14]. However, various 
cold nuclear effects such as initial parton distribution functions, Cronin enhancements and 
cold nuclear energy loss etc, are also shown to be not negligible 15|, ll6l]. Their influences 
depend on pt (see Figure [I]), i.e. shadowing effects for low px suppression, Cronin effects 
for intermediate px enhancement, isospin and EMC effects for high px suppression. There is 
also some room for a cold nuclear matter energy loss, which can reduce the overall yield [l3] . 



Dihadron correlation techniques have certain advantages over single particle observables 
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FIG. 1: Nuclear modification factor, RdAu, for charged hadrons from RUNS d+Au UM. Various 
cold nuclear effects responsible for the deviation from unity at different pT are indicated. 



in constraining cold nuclear effects. While single hadron spectra contain all contributions 
(including low Q^, non-perturbative processes), dihadron correlations are more sensitive to 
hard-scattering processes. Previously published RUNS d+Au results indicate little modi- 
fication of jet properties [l^. However, those results are focused primarily on a high 
region that might not be very sensitive to cold nuclear effects. In this manuscript, we ex- 
tend the measurements to low and intermediate pt, where the cold nuclear effects are more 
pronounced. This study allows us to have a better understanding of the final state effects 
in Au+Au collisions. 

Pairs in correlation analyses at RHIC energies are usually decomposed into a jet part and 
underlying event (UE) part (see Figure [2]). The jet part contains correlated pairs from jet 
fragmentation and associated medium response, the UE part includes the combinatoric pairs 
that are uncorrelated with the trigger particle. Generally speaking, a rigorous decomposition 
is already impossible in p+p collisions, due to non-perturbative, long range correlations 
intrinsic to QCD. In A+A collisions, the situation is further complicated by strong coupling 
between jet and the flowing medium. The decomposition only makes sense if the disturbance 
caused by the jet is localized; but if the disturbance is dissipated to the whole medium, then 
all particles are correlated with each other and the decomposition becomes very difficult and 
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FIG. 2: Schematic illustration of dihadron azimuthal distribution in for p + p and d + Au collisions. 
It has two peaks corresponding to near- and away-side jet, and a component representing the 
underlying event pairs. 



highly model dependent. This problem is not unique to dihadron correlations, it is also a 
serious issue for multi-particle correlations and analyses requiring full jet reconstruction. 

Because of the intrinsic connection between jet and the UE, it is important to study both 
of them simultaneously in heavy-ion collisions, and to understand their mutual influences. 
Admittedly, this is a difficult task. As a first step in addressing this problem, we investigate 
the UE in simpler systems, i.e. p+p and d+Au, which do not require elliptic flow subtraction. 
Our studies not only provide the necessary handles on the cold nuclear effects, but also fulfil 
longstanding interests of high energy physics community 0, 21 1. In p+p collisions, the 



UE is one of the most important backgrounds for QCD processes and for the Higgs boson 
search 221]. The study of UE at RHIC energy can provide important inputs for tuning the 
y/s dependence of phenomenological models, such as PYTHIA |23| and HERWIG |24l |. 



II. ANALYSIS 

This analysis uses the minimum bias p+p data from 2005 (2 billion events) and d+Au data 
from 2008 run (1.6 billion events) at a/snn = 200 GeV. A standard event mixing technique 
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is applied to correct for finite pair acceptance in azimutli, and a ZYAM procedure f^] is used 
to decompose the jet function into a jet part and a A0 independent pedestal. Since the jet 
signal to background ratio is large in most kinematic region considered here and jet width 
is rather narrow (except at very low px) , this should be a rather safe method for estimating 
the jet yield. 

In most jet correlation analyses, it is customary to use the per-trigger yield, PTY = 
TMggerYieid ' mcasure the jet multiplicity. To quantify jet modification, we usually compare 
jet yield with p+p collisions using one of the three nuclear modification factors: 

1. Modification of trigger yield 

TriggerYielddAu 



RdAuiPr) 



Ncou X Trigger Yieldpp 
2. Modification of pair yield 

JdAuiPT^Pr) 



^ 6 , PairYielddAu 



Ncou X PairYieldpp 
3. Modification of per-trigger yield 

dAu KPT, Pt) Per-Trigger Yieldpp 

where we use superscript "a" and "b" to indicate the two particles in the pair. In the absence 
of nuclear effects, both pair yield and trigger yield should scale with Ncou, hence deviations 
of JdAu and RdAu from unity can be attributed to nuclear effects. 

It is straightforward to show that the three quantities are related to each other via the 
following relation: 

JdAu(pT;PT) = IdAu (Pt5 Pt) R-dAu (Pt) (l) 
= IdAu (PT'Pt) R-dAu (Pt) 

This equation tell us that I^au 7^ 1 niay not necessary imply that jet is modified, it could 
be due to modification of trigger yield. This could happen if some trigger particles do come 
from jets and the yield of these triggers are modified. In this case, the jet multiplicity is 
simply re-scaled by a constant factor relative to p+p collisions, but the jet shape should 
remain unchanged. For this reason, IdAu may not be a good variable for low px correlation. 
Instead, we should use JdAm which directly reflects the jet yield modification. 
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FIG. 3: Per-trigger yield distributions for fixed pT bins from low pT to high pT- It demonstrates 
the disappearance of the modifications towards higher px- 



III. JET MODIFICATIONS IN D+AU COLLISIONS 

Modification of jet properties can be reflected by both the jet shape and jet yield. Figure[3] 
shows a sample of per-trigger yield distributions for p+p and d+Au collisions. A clear 
enhancement of the amplitude is observed at both the near-side and the away-side for d+Au 
collisions. The increase is limited to the low pt region and disappears when trigger and 
partner px values rise above 2-3 GeV/c. This enhancement of per-trigger yield is about 
factor of two at low px as measured by IdAu in Figure H]). It gradually disappears at large 
trigger and partner px- 

As we argued before, a fraction of the enhancement seen in Fig. H] is due to the suppres- 
sion of trigger yield at low px (see Fig. [T|, which increases the per-trigger yield. To avoid 
that, we measure instead the absolute jet pair yield per event and construct the pair yield 
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FIG. 4: IdAu for the near-side and the away-side for 0-20% centrality selection. 

nuclear modification factor, JdAu, for each combination and p^. A compilation of JdAu 
measurements are shown in Figure [5] for central d+Au collisions. It is plotted as a function 
of pair proxy energy p^™ = Pt + Pt' separately for the near- and the away-side. The high 
Pt — h correlation data from RUNS d-|-Au collisions [3| are included in this compilation. 
Figure [5] shows that the JdAu values approximately follow a common curve. It increases with 
p^™ at low pt, and peaks at a level significantly above one around p^™ ~4 GeV/c, then 
decreases towards larger p^™. We do not see much yield modification for peripheral d-|-Au 
collisions (Figure [6]). The shape of the enhancement of Figure [5] resembles the Cronin-like 
peak seen in single particle production. However, the level of enhancement is much bigger 
and it does not exhibit shadowing-like suppression at low pt- 

To study the modification of jet shape, we fit the jet yield distribution with a double 
Gaussian function and extract the near- and away-side width. Some examples of such fit are 
shown in Figure [31 The summary of the Gaussian widths from the fit are shown in Figure [71 
The near-side widths are identical between p-|-p and d-l-Au. The away-side widths indicate 
a small broadening in central d+Au collisions, which is expected from multiple scattering 
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FIG. 5: JdAu in 0-20% d+Au centrality selection versus p^™ = Pt + Pt near-side (top panel) 
and away-side (bottom panel). 



effects. But the level of broadening is well within the quoted 15% systematic errors. Work is 
underway to refine the systematic errors, such that we can better quantify this broadening 
in the future. 
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FIG. 7: Near-side width (left panel) and away-side width (right panel) obtained from double 
Gaussian fit in central d+Au (0-20%), peripheral d+Au( 60-88%) and p + p collisions. 



IV. THE UNDERLYING EVENT IN P+P AND D+AU COLLISIONS 

Having examined the jet shape and yield modifications, we proceed to study the properties 
of the underlying event (UE). A natural question to ask is what is the relative contribution 
from the jet and the UE. The answer to this question gives us a first order estimation 
of the contribution from hard-scattering processes in p+p and d+Au collisions. And it 
is essential for us to properly understand the modification of these processes in Au+Au 
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FIG. 8: Fraction of pairs in p+p collisions contained in the jet peak, integrated over the full 27r 
range. The systematic errors shown reflect the ZYAM uncertainties. 



collisions. Figure [8] shows the fraction of pairs in the jet peak (jet pair fraction) as a 
function of pt in p+p collisions. The jet pair fraction is already more than 15% at lowest pt 
bin (py ^ 2 and ~ 0.5 GeV/c) and quickly increases to close to 100% at high pt- This is 
much bigger than the typical level (few percent) seen in central Au+Au collisions. We want 
to point out that this estimation only gives the lower limit of jet contribution, since the UE 
contains contributions from uncorrelated jets as well as those from large angle radiations. 

In heavy ion collisions, the jet signal is difficult to extract because of the large UE level. 
In p+p and d+Au collisions, we face the opposite problem: jet signal is so big that the 
UE can be strongly influenced by effects associated with the hard-scattering process such 
as the initial and final state radiation. Our definition of the UE is different from previous 
approaches used by the CDF Collaboration [19], where the azimuth correlation is made 
between reconstructed jets and charged hadrons. This method is cleaner in separating the 
jet from the UE. However, the dihadron correlation method is still useful for two reasons. 
First, full jet reconstruction is problematic and questionable at px < 10 GeV/c, where 
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dihadron correlation method can still be used. Second, it is currently the only method 
which allows a systematic study in p+p, d+Au and Au+Au at RHIC. In order to avoid 
potential confusion, we thereafter refer the UE obtained by dihadron correlation method as 
"pedestal". 

Due to the need for event mixing to account for finite azimuthal acceptance, PHENIX 
usually defines the correlation function as the ratio of pair distributions from same-events 
and mixed-events 9|, where each is normalized separately by the number of events: 

UEpairs fr\\ 

In the second part of the equation, we decompose the foreground into the jet part and 

pedestal part via the ZYAM approach, and we use the fact that mixed event yield equal to 

the product of trigger yield, (n^), and partner yield, {rih). Because this way of constructing 

the correlation function, it is convenient for PHENIX to measure the pedestal yield relative 

to the p+p single particle yield. We define a ratio (: 

^ UEpairs/ (na) Assoc. Pcdcstal Yield Per-trig . 
{uh) Min. Bias yield Per-event 

The advantage of this quantity is that the tracking efficiency cancels in the ratio. One can 
simply multiply ( with published single particle yield 2^ to obtain the pedestal yield. 

Figure [9] summarizes ( values for various trigger pt bins as a function of partner pt- The 
( values are always above one and increase strongly with the p^ of the two hadrons. The 
change with p^ may be due to increase of the initial and final state radiation. The fact that 
> 1 can be attributed to a centrality bias in p+p events, caused by the requirement of a 
high Pt hadron pair. The impact parameter for these events are usually smaller than that 
for the minimum bias events used for the mixing, which may lead to more underlying event 
activity due to multiple parton-parton interactions (MPI). 

Figure [TO] compares the pedestal yield per-trigger integrated from 0.6 < p^ < 5 GeV/c 
as a function of trigger pt- The pedestal yield is expressed as transverse density per unit of 
azimuth and pseudo-rapidity, l/NadN/d{Ai]A(f)), similar to [13]. The integrated UE yield 
increases rapidly at low trigger px, then increases more slowly at high trigger px- Again, 
the initial increase is related to the increase of MPI, while the slower rise for larger trigger 
Pt is mainly due to the initial and final state radiation effects. 

Because of its non-perturbative nature, our current understanding of the UE is mainly 
obtained through tuning phenomenological monte carlo event generators, such as PYTHIA 
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FIG. 9: The summary of the C, in various trigger pT bins as function of associated pT in p+p 
cohisions. 



and JIMMY 26|], to match the data. In this analysis, we compare our results with calcu- 
lations from PYTHIA. PYTHIA has a large set of parameters that can affect the jet yield 
and the UE yield. To facilitate the usage, PYTHIA pre-packages several collections of de- 
fault parameter values, known as PYTHIA tunes (since version 6.410), each is identified 
with a unique integer number |27|]. One of the popular tunes is the Rick Field tune that 
can reproduce the CDF RUN2 data, known as TUNE A (100). Another popular tune is 
TUNE SOA, which is based on the new PYTHIA UE framework introduced since version 6.3. 
The comparison also depends on the jet fragmentation scheme because we do the azimuthal 
correlation at the hadron level instead of reconstructing the full jet. In this analysis, we 
considered the following three different settings from PYTHIA 6.419. 

• TUNE A (100), with parameter set tuned to CDF by Rick Field, string fragmentation. 

• TUNE A (100), but with independent fragmentation. 
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FIG. 10: Pedestal yield per-trigger per unit of azimuth and pseudo-rapidity as a function of trigger 
Pt- The shaded band indicates the yield from minimum bias p+p yield (corresponds to 21.9mb). 



• TUNE S0A(303), with New UE/MI framework, string fragmentation. 

For each setting, we generate enough simulated events, repeat the same ZYAM procedure 
as for the real data analysis, extract the jet yield and pedestal yield and compare with our 
measurements. A satisfying parameter setting should be able to reproduce all observables : 
single hadron spectra, jet yield and pedestal. Figure [TIlfT3l show such comparison of the data 
to all three settings. None of the settings can describe all three observables simultaneously. 
In general, TUNE A with string fragmentation over-predicts the jet yield, but has the best 
match for the pedestal yield; TUNE A with independent fragmentation does a better job 
for the jet yield but it does a equally poor job for single spectra; TUNE SOA with string 
fragmentation can describe the single spectra and jet yield reasonably well, but greatly 
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under-pr edicts the pedestal yield ^. These comparisons show that our data can be used to 
optimize the parameter values at RHIC energy, and to tune their dependencies on collision 



^ We would like to point out our UE study focuses at low px {< 5 GeV/c), it has been shown that the UE 
at high pt in general has better agreement with the PYTHIA tunes [28 1 
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FIG. 13: Comparison of the C values with PYTHIA tunes in three trigger px ranges. 

energy. 

If the enhancement of the pedestal is due to MPI (whose probability is proportional to 
overlap function), we should expect a strong centrality dependence of the pedestal yield in 
d+Au collisions. A systematic study of the pedestal yield from p+p to d+Au collisions thus 
can help us to understand the nature of the MPI in p+p. Figure [TH shows the pr dependence 
of ( values for p+p and d+Au collisions. The ( value at fixed px decreases towards central 
d+Au collisions. The central d+Au collisions also show much weaker increase of ( value 
with px- These results suggest that the centrality bias due to the triggering condition is 
much smaller in d+Au collisions than in p+p. In this case, the pedestal yield should be very 
close to the single particle yield for the corresponding centrality bin. 

A simple model is used to understand the scaling behavior of the pedestal yield. In 
this model, we assume that hard-scattering giving rise to the jet signal occurs only in one 
nucleon-nucleon collision in each foreground d+Au event, all other nucleon-nucleon collisions 
in the same event are assumed to be the same as the normal minimum bias p+p collision. 
In this case, the ambient particle production should scale as RdAu- The pedestal yields in p 
+ p and d + Au, UEjau and UEpp, should be related to each other through the following 
relation. 

UEdAu = UEpp + RdAu(Ncoii - l)Min.BiasYieldpp (4) 

From this equation, we can derive the following relation relating the ( values in d+Au and 
p+p 

CdAuNcoll = Cpp/ Cmb + RdAu(Ncoll — 1) (5) 
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FIG. 14: The Centrality dependence of the (" in d+Au colhsion and compared with p + p. The 
dashed Unes are the estimated pedestal level based on Eq. [5l From top to bottom, the fom' lines 
corresponds to 60-88%, 40-60%, 20-40% and 0-20% d+Au colhsions, respectively 



where emb=0.694 is the correction factor accounting for p+p trigger efficiency in PHENIX. 
The calculated ( values for d+Au collisions from this formula are indicated by the curves 
in Figure [TH Our model under-predicts the data at high and in central d+Au collisions, 
suggesting that more activities for ambient nucleon-nucleon collisions are preferred. It may 
also imply that initial state radiation associate with hard-scattering can further scatter with 
the Au nuclei, thus further increase the pedestal yield. 



V. DISCUSSION 



The results shown in Figure [MZI are intriguing: while the jet shape in d+Au collisions is 
not modified much, the yield of correlated pairs is significantly enhanced relative to binary 
scaled p+p collisions. This observation seems to imply that hard-scattering cross section is 
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enhanced but jet fragmentation and dijet acoplanarity are not modified much. Furthermore, 
the underlying event level in d+Au collisions exceeds their corresponding single particle 
yield (Figur dT^ . This means that the cold nuclear matter effects already contribute a 
significant fraction of yield enhancement seen in Au+Au, but the shape modification in 
Au+Au collisions is mostly due to final state effects in dense medium. To make sense of 
these results, we examine in the following the list of known effects in d+Au collisions, and 
speculate the roles they may play in dihadron correlation results. 

• Shadowing/CGC effect: This effect is responsible for the suppression of the single 
particle yield at low px- But our data seems to imply that it does not suppress the 
hard-scattering processes in the kinematic range considered in the analysis. 



Power corrections 29|: The expectation of binary scaling relies on the factorization 
theorem, which should break down at small . The lowest pair momentum in our 
analysis is at around 1-1.5 GeV/c, which corresponds to a value starting at around 
~ 1 — 2(GeV/c)^. At such small value, the modification jet signal due to power 
corrections might be significant. 



Multi-parton collisions [30|, IMI : This refer to the situation where more than one hard- 
scattering happens in a given event, such processes in general break the factoriza- 
tion. The simpliest case is double-parton collisions (two independent hard-scattering 
occurs in same event). The rate such scattering is greatly enhanced in nuclear envi- 
ronments. A simple estimation show that ratio of double-parton scattering to single 



hard-scattering is cxf /(xf ^ 0.5(A/10)°-^ [3 

Initial state multiple scattering effects: The rate of multiple scattering, where a hard 
parton undergo one hard-scattering plus several soft scattering, is significantly en- 
hanced in p+A collisions. Combined with steeply falling spectra can shift initial 
parton to higher px-, it can increase the jet pair yield and pedestal yield. The fact that 
near-side jet shape is not modified suggests that the multiple scattering happens at 
the parton level before fragmentation. 

Interaction of initial state radiation with the remaining nuclei: In p+p collisions, such 
radiation simply fragments into final state hadrons. But in the d+Au environment. 
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they may undergo additional scattering with the nuclei, thus further increasing the 
observed jet pair yield and pedestal yield. 

• Cold nuclear jet energy loss: This effect can increase the yield of the soft partons, 
which in turn leads to more pair yield at low pt- However, this scenario also leads to 
a suppression of the high pt pair yield, which is not seen in our data. Thus this effect 
alone can not explain our data. 

• Isospin and EMC effects: It can affect jet pair yield at very high pj^"^ where the quark 
jet contribution are important. In the kinematic region of our analysis, most jet pairs 
should be dominated by gluon jets, hence these effects should not be important. 

The exact contributions from each of these effects are not clear. It is possible that several 
effects conspire to give the observed modification on the jet correlation and the UE. Future 
theoretical efforts should use both single spectra and dihadron correlation results to better 
constraint these effects. 

In summary, we performed a detailed study of the dihadron correlation in d+Au collisions 
over a broad pt range. A large Cronin-like enhancement is seen at low pair proxy energy 
while the jet shapes show little difference from p+p. This is the same region where a large 
modification of the jet properties are also observed in Au+Au coUisions. This observation 
favors initial state multiple scattering effects at the partonic level. The jet properties in p+p 
are compared to several PYTHIA calculations. The most popular PYTHIA tunes that work 
well at Tevtron energy fail to simultaneously describe the single particle yield, jet pair yield, 
and especially the pedestal yield at RHIC energy. The scahng behavior of the pedestal yield 
from p+p to d+Au suggests that underlying event yield in d+Au exceeds a simple sum of 
one hard-scattering p+p collision and Ncou — 1 minimum bias p+p collisions. Theoretical 
investigation of the centrality dependence of the underlying event in d+Au can help us to 
understand the nature of the multiple parton interaction in p+p and d+Au. 
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